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Abstract 

We revisit the idea of generating the Higgs /i parameter through a spontaneously broken 
Peccei-Quinn (PQ) symmetry in gauge-mediated supersymmetry breaking scenario. For the 
messenger scale of gauge mediation higher than the PQ scale, the setup naturally generates 
fj, ~ wisoft and the Higgs soft parameter B < (9(msoft) with the CP phase of B aligned to the 
phase of gaugino masses, while giving the PQ scale vpQ ^ -y/^^softA, where m^oit denotes the 
gauge-mediated gaugino or sfermion masses and A is the cutoff scale which can be identified as 
the Planck scale or the GUT scale. The PQ sector of the model results in distinctive cosmology 
including a late thermal inflation. We discuss the issue of dark matter and baryogenesis in 
the resulting thermal inflation scenario, and find that a right amount of gravitino dark matter 
can be produced together with a successful Affleck-Dine leptogenesis, when the gravitino mass 
7713/2 = C»(100) keV. 
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I. INTRODUCTION 



Weak scale supersymmetry (SUSY) is one of the most attractive candidates for new physics 
beyond the standard model (SM) at the TeV scale It provides an appealing solution to 
the gauge hierarchy problem, and also the successful unification of gauge couphngs at the scale 
Mqut ~ 2 X 10^^ GeV. On the other hand, the absence of unacceptably large flavor or CP 
violations requires a rather special type of supersymmetry breaking which yields flavor and CP 
conserving soft terms. Supersymmetry breaking through gauge mediation ^ provides flavor 
conserving (possibly CP conserving also) soft terms in a natural manner as the structure of 
soft terms is determined mostly by the SM gauge interactions. One potential difficulty of 
gauge mediation mechanism is the generation of the Higgs fi and B parameters having a right 
size for the electroweak symmetry breaking. If the Higgs sector communicates with the SUSY 
breaking sector to generate fi ~ "n^soft, where msoft denotes the gaugino and sfermion masses 
in gauge mediation, one often finds B ~ 87r^msoft, which is too large to achieve a successful 
electroweak symmetry breaking. There have been many attempts to solve the n problem in 
gauge mediation, including those in Ref. j^. 



As was noticed in the original work of Kim and Nilles |4|, a satisfactory solution of the 
problem should provide a theoretical reasoning for the absence of the bare fi term with /i ~ A, 
as well as a dynamical mechanism to generate fi ~ rrisoft together with B ~ msoft at the weak 
scale, where A denotes the cutoff scale of the model which can be taken as either the reduced 
Planck scale Mp, ~ 2 x 10^^ GeV or the GUT scale Mqut ~ 2 x 10^^ GeV. It was further 
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noticed in Ref. that the Peccei-Quinn (PQ) symmetry solving the strong CP problem might 
play a crucial role for the fi problem as well. The U{1)pq symmetry might forbid the bare /i 
term, while allowing the following non-renormalizable term in the superpotential 

—X'^HuHd, (1) 

where X is a PQ charged SM singlet field whose vacuum value breaks U{1)pq spontaneously. 
If the PQ sector of the model couples to the SUSY breaking sector to stabilize X at 



VPQ 



(X) ~ VmsoftA, (2) 



and the F-component of the stabilized X satisfies /X < C(msoft), the resulting /i and B 
(at the weak scale) have a right size for successful electroweak symmetry breaking. 



2 



The Kim-Nilles mechanism was discussed originally in the context of gravity mediation with 
Later it was realized that the mechanism can be implemented also in gauge 
mediation [r]. However, the specific models discussed in \}\ involve U{1)pq which is assumed 
to be an i?-symmetry. In such models, the nonzero vacuum value of the superpotential, which 
is required to tune the cosmological constant vanish, should appear as a consequence of the 
spontaneous breakdown of U{1)pq, and this makes a complete realization of the setup quite 
complicate. In this paper, we revisit the Kim-Nilles mechanism to generate fi and B in gauge 
mediation, while focusing on the case that U{1)pq is not an i?-symmetry, but an ordinary 
anomalous global symmetry. It is noticed that such class of models can have a distinctive 
cosmological feature such as a late thermal infiation triggered by the PQ sector jsj. We then 
need a late baryogenesis after thermal infiation as well as a mechanism to produce a right 
amount of dark matter. We find that a right amount of gravitino dark matter can be produced 
after thermal infiation when vpq = C(10^ - 10^°) GeV and mg/s = C(IOO) keV. With the 
nonrenormalizable term ([T]) and also the seesaw term for light Majorana neutrino masses, the 
model can accommodate also a successful Affleck-Dine leptogenesis proposed in 

This paper is organized as follow. In section 2, we discuss the Kim-Nilles mechanism to 
generate fi ~ msoft and B < (9(msoft) with a spontaneously broken U{1)pq in gauge mediation. 
Section 3 discusses the cosmological aspects of the model, including the mechanisms to generate 
the right amount of dark matter and baryon asymmetry, and the conclusion will be given in 
section 4. 



II. THE KIM-NILLES MECHANISM IN GAUGE MEDIATION 

In (minimal) gauge mediation scenario, SUSY breaking is mediated by SM gauge-charged 
messengers $, ^'^ which couple to SUSY breaking field Z = M + FO"^ in the superpotential. 
Then sfermions and gauginos in the minimal supersymmetric standard model (MSSM) get soft 
SUSY breaking masses 

through the loops involving the messenger fields $, In order to implement the Kim-Nilles 
mechanism to generate the term, we introduce additional SM singlet but PQ charged super- 
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fields which break U{1)pq spontaneously, and also extra vector-like quark superfieldi*] which 
have a Yukawa coupling to the t/(l)pQ-breaking fields. With U{1)pq, one can forbid renormal- 
izable superpotential term of the t/(l)pQ-breaking fields, while allowing a nonrenormalizable 
term suppressed by the cutoff scale A of the model. Then due to the SUSY breaking ef- 
fects mediated through the Yukawa coupling to extra quark superfields, the ?7(l)pQ-breaking 
fields are destabilized at the origin. On the other hand, the supersymmetric scalar potential 
originating from the nonrenormalizable superpotential prevents the runaway of the U{1)pq- 
breaking fields, and stabilize them at an intermediate scale vpq ^ -y/^softA. With A presumed 
to be the GUT scale or the Planck scale, this scenario naturally generates a QCD axion scale 
vpQ = (9(10^ — 10^^) GeV, as well as a correct size of /i ~ Vpg/A ~ frisoit via the Kim-Nilles 
mechanism. Furthermore, in this setup one can easily obtain the Higgs B parameter at the 
weak scale which is (at most) comparable to msoft and has a CP phase aligned to the phase of 
gaugino masses. 

As a specific model to realize this scenario, we consider the superpotential 

+XzZ^^'' (4) 

where the first line denotes the usual Yukawa couplings between the Higgs fields and the quarks 
and/or leptons, including the term for small neutrino masses which might be generated by the 
seesaw mechanism [llj] with a right handed neutrino mass M^v far above the weak scale. Here 
the fiavor indices are omitted, and i/d, He and yi, should be understood as 3 x 3 matrices. 
As we will see, the above model has a variety of interesting cosmological features, including a 
late thermal infiation associated with the PQ phase transition in the early universe. Although 
we do not specify the origin of M^v here, an interesting possibility is that is generated as 
a consequence of U{1)pq breaking, which would give M^r ~ {X), so that the seesaw scale is 
identified as the PQ scalqll 12(]. As we will see in the next section, such setup can be useful 
also for a successful Affleck-Dine leptogenesis after thermal infiation. 

The second line of the superpotential (jl]) is the PQ sector generating the Higgs /x parame- 
ter through the Kim-Nilles mechanism, while providing a QCD axion to solve the strong CP 
* To keep the successful unification of gauge couplings in the MSSM, these extra vector-like quarks can be 

extended to form a full GUT multiplet. 
^ Of course, then the Yukawa couplings between and the left and right handed neutrinos should have 

appropriately small values to produce the observed^eutrino mass-square differences and mixing angles. 



problem 13|]. The third hne is for the minimal gauge mediation of SUSY breaking, where Z is 
the SUSY breaking field with 



(Z) = M + F9^. 



(5) 



Note that we can always make Ax, ^i, K2 and \z{Z) all real and positive through appropriate 
field redefinitions, and we will take such field basis in the following discussion. To be specific, 
we also assume that the cutoff scale A is around the GUT scale 



A ~ M, 



GUT 



2 X 10^^ GeV. 



(6) 



Although different choice of A would change the value of the PQ scale, the fi parameter obtained 
by the Kim-Nilles mechanism is independent of A and always of the order of rrisoft as long as the 
dimensionless parameters ni and K2 have a similar size. Note that the superpotential @ takes 
the most general form (up to dim = 4 terms) allowed by the SM gauge symmetries, i?-parity 
and U{l)pQ, where the U{1)pq charges are given as follows: 



Field 


Z 


X 


Y 


Hu 












<|)<|)C 


PQ charge 





1 


-3 


-1 


-1 


1 


1 


1 


-1 






Let us now discuss the vacuum configuration of the f/(l)pQ-breaking fields X and Y. As 
the gravitino mass is much smaller than the weak scale, we can safely ignore the supergravity 
effects. Then the scalar potential of X, Y can be well approximated by the global SUSY 
potential including the soft SUSY breaking terms induced by radiative corrections: 



V{X,Y) = V,,i,{X) + 



"^1 |X|6 + ^|X|^|FP 



36A2' 



4A2' 



(7) 



where 



l^soft(X, Y) = m\\X\^ + ml\Y\^ + {a^^^X^Y + h.c.) . 

If the messenger scale M$ of gauge mediation is above the PQ threshold scale 

M^ = \z{Z)>\x{X), (8) 

which is in fact necessary to generate /x ~ msoft independently of the value of A, the soft 
mass mx at scales below M$ is generated mostly by the renormalization group (RG) running 
triggered by the Yukawa coupling XxX^^'^. The RG equation for m\ is given by 



(iln/i^ 



3Ai 
167r2 



(9) 



where the factor 3 is the color factor, and m^,m^c and v4x*^<= are the gauge-mediated soft 
scalar masses and trilinear scalar coupling, respectively, for the squark components of \E', \E'^. 
As mx and at the messenger scale are negligible compared to niq, ~ msoft, the soft mass 

mx at lower renormalization point {X) is determined as 

m|(|X|)^-K+m|j|^ln|^, (10) 

where we have ignored higher powers of g^ln(M$/|X|) with the assumption that M$ is not 
so far above (X). Since the PQ breaking scale (X) is constrained to be of O{10^ — 10^^) GeV, 
while the messenger scale should be lower than O{10^^) GeV in order for the gauge mediation 
to give dominant contribution to soft terms, the value of ln(M$/|X|) can not be so large, and 
therefore our assumption is justified. 

In our approximation, and m^c in ( ITOj) can be regarded as the soft squark masses at the 
messenger scale M$, which are given by 



2 2 OiV<J) 



~ "^soft, (11) 



167r2 M 

where iV$ is the number of messenger pairs in the fundamental representation. Minimizing the 
scalar potential ([7]) with the tachyonic im?x given by ( fTOj) . we find 



2 /|^|2\ 3AxVln(M$/t;pQ) 
VpQ = {\X\ ) ~ ^ m^A = C(msoftA) (12) 



and the resulting Higgs /x parameter 



^^ ^-^^WHM./.., _l OM, (13) 

where we assumed that Ax, i^i and k,2 are all of order unity for the order of magnitude estimate 
in the last step. Note that fi is independent of the precise value of the cutoff scale A, while the 
PQ scale has a mild dependence on A. 

The VEV of X in f[T2|) generates an effective mass of Y through the term oc |X|^|y p in the 
scalar potential: 

^,W)n = 3\mxm\ (14) 

It also generates an effective tadpole of Y through the A term oc X^Y in the scalar potential, 
which is generated by the RG evolution of the wavefunction factor of X. We then find 

A,,(|x|)_ svmx'x f^^M^\\ ^^5) 



4V2n^ 8tt^ \ \X 



and therefore 



1^7^^) • (16) 



3v^|mx| 2v^ 87r3 V (1^1) 
With the above resuhs, one can compute the Higgs i?-parameter around the messenger scale 
which is given by 

^ +^)-0, (17) 



X J \\ X* 

upon ignoring the gravity mediated contribution of (9(m3/2). Note that the equation of motion 
of Y leads to the cancellation between the two contributions to B{Mip), making B{Mip) even 
smaller than 0{Ai^-^). As the B parameter at M<j, is negligible compared to mgoft, its low energy 
value is determined by the RG running from M$ to the weak scale. In our case, the messenger 
scale M$ is required to be higher than the PQ scale vpg ~ -y/'^soft ^gc/t = 0(10^ — 10^°) GeV. 
As a result, a sizable value of B can be induced at the weak scale, giving tan /3 = 10 ~ 20, and 
furthermore its CP phase is automatically aligned to the phase of gaugino masses. 

With the PQ sector stabilized as above, we can identify the mass eigenstates of the PQ 
sector fields and compute their mass eigenvalues. First of all, the PQ sector provides a QCD 
axion having a decay constant vpg and thus a mass ~ fTr^ir/vpQ, which corresponds mostly 
to the phase degree of freedom of X. It contains also three real scalars with a mass comparable 
to msoft, i-e. the saxion x which is mostly the modulus of X and two others from Y = yi + iy2, 
and two Majorana fermions fij {i = 1, 2) which form approximately a Dirac axino d = {di, 0,2). 
It is then straightforward to find 

m^~2|mx|, mj^j ~ ~ ma ~ Vslmxl, (18) 

where mx is given by (fTOj) . 

In summary, in our model the messenger scale of gauge mediation is assumed to be higher 
than the PQ scale, and then the U{1) pg-breaking fields X and Y are destabilized from the origin 
due to the tachyonic soft mass of X and the scalar A term associated with the nonrenormalizable 
superpotential term kiX^Y/6A. These soft SUSY breaking terms of X and Y are induced by 
the combined effects of gauge mediated SUSY breaking and the Yukawa coupling XxX'^'^'^. 
The supersymmetric scalar potential from the nonrenormalizable superpotential prevents the 
runaway of X and Y, and stabilizes them as 

1/2 / X 1/4 



<i^i>-i^('"(M) '''' 

where the messenger scale M$ > vpg = and A is the cutoff scale of the model. If we 

assume that Ax, i^i and K2 are all of order unity and A ~ Mgut = 2 x 10^^ GeV, while M$ is 
not so far above fpg, the mass scales of the model are estimated as 

/i ~ m,oft ~ mpQ = O(102 - 10^) GeV, 

vpQ ~ VmsoftMcuT = 0(10' - 10^°) GeV, (20) 

where mpQ stands for the masses of the PQ sector fields (other than the QCD axion), including 
the saxion and axino masses. The B parameter at the messenger scale is negligible compared to 
msoft, and therefore its weak scale value is determined by the RG evolution below the messenger 
scale, making its CP phase automatically aligned to the phase of gaugino masses. 



III. COSMOLOGY OF THE MODEL 



The model described in Sec. 2 has a variety of interesting cosmological implications. Because 
the PQ preserving field configuration X = F = is a local minimum of the effective potential 
at high temperature T ^ w^soft, it is a quite plausible possibility that X is settled down at 
the origin after the primordial inflation. Then the early Universe experiences a late thermal 
inflation jsl, 1^ before the PQ phase transition occurs, which might be useful to eliminate (or 
dilute) potentially dangerous cosmological rehcs such as light moduli or gravitinosH Since 
this thermal inflation will erase out any primordial baryon asymmetry, we need a baryogenesis 
mechanism operating after thermal inflation is over. As for the dark matter in our model, 
one can consider two possible candidates, QCD axion with a mass m^, ~ /ttJ^tt/'^pq and light 
gravitino with a mass 7713/2 ~ F/Mpi. However the PQ scale of our model is determined as 
vpQ ~ y/msohMcuT ~ O{10^ — 10^°) GeV, which might be too low to give a QCD axion 
constituting the major fraction of the observed dark matter pj,]. This leads us to focus on the 
possibility of gravitino dark matter. 

In the following, we briefly discuss the cosmological features of our model, while leaving 
more complete discussions for a separate paper |15:]. As we will see, for the gravitino mass 



^ Note that a UV completion of the model within the framework of supergravity or string theory might contain 
cosmologically harmful light moduli causing the so-called moduli problem. 
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range 

m3/2 = C(IOO) keV, (21) 

a right amount of gravitino dark matter can be produced after thermal inflation, together with 

a successful Affleck-Dine leptogenesis. As the gravitino mass is given by 

F IQtt^ M , , 

Mpi Mpi 

for the SUSY breaking spurion Z = M + FO"^, this range of 7713/2 suggests that the messenger 
scale of gauge mediation, i.e. M$ = A^M, which is presumed to be higher than the PQ scale 
vpQ, should be somewhat close to vpQ ~ 10^ — 10^° GeV. 

A. Thermal inflation 

Near the origin, the finite-temperature effective potential of the flat direction |X| is given 

by 

V{X) = Vo+ {PIT' - \m\m) |X|2 + . . . , (23) 

where Vq = 0{m?xv1,Q) is the potential energy at the origin, which is set to make the cosmolog- 
ical constant at true vacuum vanish, and Px comes from the couplings to thermal bath. Once 
the Universe were in a radiation dominated period with T > Vq^'^ after the primordial inflation 
is over, thermal inflation begins at the temperature 

Tb ~ Vq'^ ~ 10*^ GeV (24) 

and ends when |X| is destabilized from the origin at the critical temperature 

= ». (25) 

Px 

Soon after the end of thermal inflation, the Universe is dominated for a while by the coherent 
oscillation of |X| around its true minimum = fpg, which eventually decays into lighter 

particles to reheat the Universe. Since the saxion mass = (9(10^ — 10'^) GeV in our model 
(see ( fTSl) ). it can decay dominantly to the light Higgs boson pair h + h* through the coupling of 
the form fi''hh*6x/vpQ, where 6x denotes the saxion fluctuation around its vacuum. Assuming 
that is heavier than 2mh for the light Higgs boson mass rrih ~ 120 GeV, the decay rate of 
IXI is estimated as 



and then we find the reheat temperature is given by 

X -1/4 



Trh = ^-g*{TRR)j rfM'pi" (27) 

where g^iT^n) ~ 100 is the number of hght degrees of freedom at T = T^ih- The total number 
of e-foldings of this thermal inflation is estimated to be about 10 and the dilution factor due to 
the entropy release in the decay of \X\ is about O{10^^). This would be large enough to remove 

3- 



for instance the gravitinos produced before thermal inflation 

Our model has two other oscillating scalar fields which are mostly yi = Re(y) and y2 = 
Im(y). Although they have a mass comparable to rrix (see (fTSll ). their energy densities are 
suppressed by = (9((7^A|-/(87r^)^) compared to that of |X|, and therefore they do 

not give a significant impact on the cosmological evolution after thermal inflation. 



B. Affleck-Dine leptogenesis 



Thermal inflation erases pre-existing baryon asymmetry. One may think that an Affleck-Dine 
(AD) baryogenesis before thermal inflation with a very large initial value of AD fleld can produce 
enough baryon asymmetry which would survive after thermal inflation 17| . However it is known 



that the formation of Q-balls makes it difficult to realize such scenario I 
model can realize the late-time AD leptogenesis proposed in Refs. 

m 



Fortunately, our 



19|423|. 



In order for the AD leptogenesis to work, the MSSM flat direction LH^ is required to have 
a nonzero value at certain stage. In our case, this initial condition can be achieved as LHu 
has a tachyonic soft mass-square —Tn\^^ in the limit /i = 0, so unstable at the origin if the 
temperature drops below its critical temperatures Tlh^ = ^lHuI (^lHu ~ V^TnLHu ^^nd X is 
still staying at the origin. This requires 

|mx(0)| 



X 



and thus 



/3i > 



6A| 

7r2 



m 



In 



VPQ 



(29) 



(30) 



where we have used mx(0) ~ imxivpg) together with the result ffTOl) for mxivpg). In our 
model, /3x receive a contribution from the exotic quark superfields giving 



A/3_ 



4 ' 



(31) 
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It turns out that it is difficult to satisfy fl5U]) only with A/3|^ for typical parameter values of our 
model. However this difficulty can be easily avoided if the field X couples to the right-handed 
neutrinos to generate the seesaw scale |l2l| . Then there will be an additional contribution to 
I3\ from the Yukawa coupling A nXNN t2J] , 



\(j2x% + SXj,y (32) 



io| 



with which the condition fl30|) can be satisfied with a reasonable value of Atv- 

Once the key condition (I5U]) for AD leptogenesis is satisfied, the AD field LH^ rolls down 
to nonzero value at the temperature T ~ ^lHu- If ^ drops further down to T^, |X| rolls 
away from the origin to generate nonzero and then LHu gets a positive mass-square due to 
the contribution from /i^. As a result, LHu rolls back to the origin with an angular motion 
generated by CP-violating terms in the scalar potential. The lepton asymmetry associated with 
the angular motion of AD field is finally converted to baryon asymmetry through the sphaleron 
process. The resulting baryon asymmetry at present is estimated as 

ub SnLTRH 

— ~ , [66) 

s 8 UxTrix 

where Ux ~ rrixVpQ is the saxion number density for coherently oscillating saxion field \X\, 
and riL is the lepton number density associated with the angular motion of the AD field. In 
fact, ul depends on many details of the full scalar potential, including the terms associated 
with the lepton number violating neutrino mass term in the superpotential. It depends for 
instance on the initial displacement of the AD field from the origin, curvature of the potential 
in angular direction, CP phase, e.t.c. Using the results of 10|, l20|, l2l|, 1231], we find that a value 
of riL = 0(10^^ — 10^^) GeV^ can be achieved under a reasonable assumption on the involved 
model parameters, and therefore the AD leptogenesis after thermal inflation can produce the 
observed value of n^/s ~ 10~^° within the uncertainties in the involved parameters. 



C. Dark matter 



In our model the lightest supersymmetric particle (LSP) is the gravitino, because all other 
supersymmetric particles including axinos have mass of order msoft, and the gravitino mass 
is much smaller than rrisoft- Thus light gravitino is the prime candidate of the dark matter. 
On the other hand, thermal inflation also dilutes pre-existing gravitino relics. After that, by 
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decay of |X|, the Universe is reheated with temperature Trh ~ ITeV as in Eq. (1271) . At this 
temperature, most of MSSM fields are thermahzed and will produce light gravitinos. We can 
divide this process into two parts. One is thermal (TH) production in which gravitinos are 
produced by scatterings and decays of the MSSM fields in thermal bath. The other is non- 
thermal (NTH) production in which gravitinos are produced by out of equilibrium decays of 
the frozen relics such as the next LSP which is the ordinary LSP (OLSP) in the MSSM sector 
or the axino. The corresponding relic density of the gravitino can be represented by 

ny,h' = njjie + fi^™/.^ ^ 2.8 X 10^ (y^) {ysl^ + ysT) (34) 

where ^3/^^"""^ = '"'3^2™ (-^) / ^C^) yield of the gravitino which is produced by (non-)thermal 
process. At Trh ~ ITeV, gravitinos from thermal production can provide a right amount of 
cold dark matter, ~ 0.1, if the mass is jlG^ 

m3/2 ~ 100 keV. (35) 

For the NTH production of the gravitinos, the contribution from the OLSP decay is small 
enough in the above range of ^3/2, but we have to pay attention to the production from the 
axino decays. Although the axino couplings to the MSSM particles are suppressed by 1/fpQ, 
they are still large enough to generate a significant axino abundance from the thermal bath, 
f the axino is stable, one needs Trh ^ ^soft in order to suppress its relic density sufficiently 
25I . In our case of Trh ~ "^soft? thermally generated axinos may produce a large number of 
gravitinos from their decays. If axinos decay only to gravitinos, the non-thermal relic density of 
gravitinos turns out to be too large. Thus, the axino decay to the gravitino must be suppressed. 
For this, let us now consider the following axino mass range: 

+ ruh < rria < - ruh (36) 

where x is the OLSP. Then the dominant production and decay channels of the axino come 
from the higgs(/;.)-higgsino(/;,)-axino(a) coupling: 

d^e^X^H^Ha = -^hhd + ■■■. (37) 
2A vpQ 

That is, axinos are produced thermally by the process h — )■ hd, and then they decay mainly 
through a — )■ hx- Denoting the decay rate of higgsino to axino as r(/i — )■ hd), one finds the 
thermal axino abundance as follows: 



_ 135C(3) T{h hd) 
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(38) 



= w7:::r >^" (39) 



where (7* ~ 200 is the relativistic degrees of freedom and H ^ 0.33y^T^/Mp is the Hubble 
parameter at the termperature T. Then, the non-thermal abundance of the gravitino is 

.NTH _ r(a 1^3/2 a ) ^^TH 

in which partial decay rates of axinos are given by r(a — 1^3/2 a) = m\ /QQTimli^Ml, [26| with 
"03/2 being the gravitino and T{a — > xh) = O'^fi'^ma/SiTVpQ. Here 6 parameterizes the OLSP 
fraction in the Higgsino component. For our estimation, we will use 6 ~ mzsw/f^ which is 
valid in the limit of the large Higgsino and small gaugino masses. Combining (!38|) and (139|1 . 
we get 

Y-r^i.,.w-u^y('-^y('-^) (40) 

V350GeV/ V ^3/2 J \ J 
which shows that the non-thermal gravitino relic density can be safely neglected. Finally let us 
remark that the axino decays well before the OLSP freezes out, that is, r(a — > xh) > H{Tf) 
for the OLSP freeze-out temperature Tf < 20 GeV for our choice of parameters. Therefore, the 
OLSPs from the axino decay are thermalized. 



IV. CONCLUSION 



In this paper, we have examined a model to generate the Higgs jj, parameter with a spon- 
taneously broken U{1)pq symmetry in gauge mediation scenario. The PQ sector of the model 
contains U{1)pq breaking fields which have a Yukawa coupling to extra quarks. The U{1)pq 
breaking fields also have a nonrenormalizable superpotential suppressed by the cutoff scale A 
which might be identified as the Planck scale or the GUT scale. For the messenger scale higher 
than the PQ scale, the U{1)pq breaking fields are destabilized at the origin due to the soft 
SUSY breaking terms induced by the combined effects of gauge mediated SUSY breaking and 
the Yukawa coupling to extra quarks. They are then stabilized by the supersymmetric scalar 
potential from nonrenormalizable superpotential at an intermediate scale fpg ~ v^msoftA, gen- 
erating /i ~ Vpg/A ~ msoft in a natural manner. The B parameter at the messenger scale is 
predicted to be negligible, and therefore B at the weak scale is determined by the RG evolution 
below the messenger scale. 

The model has a variety of interesting cosmological features associated with the PQ phase 
transition. In particular, a late thermal inflation is a natural possibility, which would require a 
late baryogenesis mechanism. We find that a successful Affleck-Dine leptogenesis after thermal 
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inflation can be implemented within the model. We also find that a right amount of gravitino 
dark matter can be produced after thermal infiation when vpq = O{10^ — 10^°) GeV and 
= 0(100) keV, for which the messenger scale of gauge mediation is required to be not far 
above vpq. More complete discussion of the cosmological aspects of the model will be presented 
elsewhere. 
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